An investigation of integrating the finite element method (FEM) with grey system theory for geotechnical problems

Numerical simulation is very important to solve geotechnical problems. However, it is difficult to obtain required comprehensive and accurate information such as parameters, boundary conditions, and etc. In this paper, a grey distributed parameter model, which integrates the finite element method (FEM) with the grey system theory, was proposed to address the issue. The analysis of grey properties on rock and soil system was performed. The equilibrium equations, geometric equations, physics equations and related differential equations were obtained, each of the equations contains grey parameters and variables. And the discretization and solution methods of the FEM with the grey variables were discussed. An example of deep-buried circular mining tunnel was applied to test the proposed model. The calculation results were compared with those of the exact solution (analytical solution) and the classical FEM, respectively, through which the rationality of the proposed model was demonstrated. For the first time, grey variables and grey parameters are defined in geotechnical numerical simulation. The expressions of basic equations with grey variables are given respectively. A grey distributed parameter model which integrates the FEM with the grey system theory is proposed to solve geotechnical problems, and the optimal solution to the proposed model is determined through calculation and comparison of an application example. The proposed numerical model with grey variables not only has the advantage of grey system theory, but also greatly improves the adaptability and application effect of the model, which contributes to the prediction and evaluation problems in geological engineering, geotechnical engineering, water conservancy engineering and civil engineering with complex structures.


Introduction
Geotechnical engineering solves problems in stability of surrounding rock in various underground chambers (such as tunnels, building foundations, goaf, and etc.) and also in the formation and change of various geological disasters on earth's surface. The empirical formula a1111111111 a1111111111 a1111111111 a1111111111 a1111111111

OPEN ACCESS
Citation: Li Y, Zhu G, Zhang Q (2022) An investigation of integrating the finite element method (FEM) with grey system theory for geotechnical problems. PLoS ONE 17(6): e0270400. https://doi.org/10.1371/journal. pone.0270400 method, engineering analogy method and numerical simulation method are used for solving geotechnical problems nowadays, among which the numerical simulation provides evaluation results of surrounding rock stability more scientifically. However, these methods are based on a 'white' (a well-defined) system. In fact, the structure of geotechnical system is complicated, and most is in the blind area of people's vision. Data of rock and soil system obtained by various engineering exploration methods are limited and partial, and it makes the geotechnical system a non-essential grey system which could be better studied with grey system theory [1].
In the field of geotechnical engineering, grey system methods have been actively applied and some achievements have been made. LU and Rosenbaum [2] predicted future ground movement based on geotechnical properties and historical behavior with artificial neural networks (ANN) and grey system theory. Jing et al. [3] established a simulation experiment based on the grey system theory with the actual measurement data of the studied surface. The vertical displacement of the surface was used as a reference array, and a comparison array with four groups from the experiment constituted a grey relative analysis system. The relationship between different mining subsidence and vertical displacement of the surface was determined, which provided a theoretical basis for the layout of the mine work area and the reduction of surface building damage. Bachir et al. [4] applied grey system to predict working face pressure, with particle swarm optimization (PSO) algorithm obtaining optimal parameters in the grey model, as a result, achieving good prediction effect. Deng et al. [5] applied the grey system theory on predicting crown settlement of Jiulianshan tunnel, and compared it with the traditional regression analysis method. Lei et al. [6] established the classification table of surrounding rock stability according to the geological data. On the basis of grey relational analysis and fuzzy recognition principle, the objective function was constructed by the least square method. The relational analysis of five indexes affecting the stability of surrounding rock was carried out, and the grey optimal classification model to evaluate surrounding rock stability was established. Tasci and Kose [7] used multivariable grey models GM (0, N) and GM (1, N) to determine the deformation on the crest of the Keban Dam in Turkey. Miao et al. [8] used the Verhulst-Grey model to calculate the parameters of the landslide time forecast model established on the basis of displacement monitoring data. Zhang et al. [9] proposed a model which combined the optimized grey discrete Verhulst model with back propagation (BP) neural network to better predict foundation pit settlement. Zhang et al. [10] updated the traditional grey model (TGM) (1, 1) to the TGM (1, 1, p, q) grey model by adding a nonlinear time correction factor and optimizing the weighting coefficient of the background value. Wang et al. [11] proposed a grey neural network model based on the grey theory with the neural network algorithm, and analyzed the roof compaction data of the return air working face in a mine in Xuzhou. It not only eliminated the shortcomings of neural network model, but also made up for the lack of self-feedback adjustment of grey network. MATLAB was used to simulate and test the roof pressure prediction results of the model, and the results were compared with those of single grey model and single BP neural network model. As can be seen from the work cited above, the application of grey system theory in geotechnical problems mostly employed grey models with lumped parameters which can only predict and evaluate the studied object as a whole, making the simulation of any point with the methods not satisfying.
Moreover, the classical numerical simulation method with distributed parameters is capable of simulating the state of any point. It was used for analysis of the stress, displacement and failure of underground engineering, slope and dam since 1960s and has made great progress with the development of computer technology. Due to the difficulty in obtaining sufficient and comprehensive data such as mechanical parameters, loads and boundary conditions, the application of the classical numerical simulation method is also limited. Therefore, studies on integrating the FEM with uncertainty methods are getting intensive. Kleiber and Hien [12] proposed a new stochastic FEM based on the perturbation method to analyze the structural response. Ghanem and Spanos [13] have proved that the spectral stochastic finite element method could forecast a variety of uncertainties in calculating system responses. Griffiths and Fenton [14] used a random FEM (RFEM) for stochastic analysis of slopes by considering a strength reduction technique. Chaudhuri and Sekhar [15] have adopted the stochastic FEM based on perturbation technique to analyze the flow and transport in a three-dimensional porous formation. Ahmed [16,17] conducted a stochastic seepage analysis through earth dams and under water-retaining structures by FEM. Zheng et al. [18] used fuzzy FEM to analyze the impact of cold wave on face slab cracking of a concrete-faced rockfill dam (CFRD). A simplified procedure based on FEM with the random field theory was developed by Hongwei Huang et al. [19], for analyzing the longitudinal performance of shield tunnels considering the longitudinal variation of geotechnical parameters. Johari and Talebi [20] considered the uncertainties of soil parameters leads to an uncertain wetting branch of the soil-water retention curve (SWRC), and made a computer program to determine the reliability indices of unsaturated slopes under different rainfall intensities and seepage flow using random FEM (RFEM). Chen and Liu [21] developed a stochastic FEM to solve the calculation precision deficiency caused by spatial variability of dam compaction quality. Liu and Wang [22] proposed a probabilistic method for simulating the entire process of rainfall-induced landslides considering spatial variability of soil properties. A two-stage method containing FEM and material point method (MPM), both with hydro-mechanical coupling, is proposed to simulate a landslide from triggering by rainfall to post-failure large deformation of soils. It can be seen from the work cited above, that FEM is integrated with many uncertainty methods, such as fuzzy-FEM, random-FEM and stochastic-FEM. However, there are few researches on integrating FEM with grey system theory, even though grey system has the advantages of solving system problems in poor information situation.
In summary, up to now, studies of the grey system theory and methods in the field of geotechnical engineering mostly employed grey models with lumped parameters, and the simulation results are lack of accuracy especially when performing analysis of any point. In terms of the classical method (analytical or the numerical) with distributed parameters, although it can be used for simulation of any point, the parameters and other information involved are treated as deterministic, making the response obtained also deterministic. But the actual geotechnical system is not completely this way, and grey system has the advantages of solving system problems with poor and uncertain information. Therefore, in this paper, a grey distributed parameter model (numerical model) which integrates the FEM with grey system theory was proposed to study geotechnical problems, and it can be used to simulate displacement of any point without deterministic input information. The grey properties of rock and soil system are analyzed; various equations with grey variables are established including mechanical equilibrium equations, geometric equations, physics equations, and related differential equations. Then, the discretization and solution method of the finite element equation with grey variables is discussed. The proposed model (FEM with grey variables) was used to analyze an example of deep-buried circular tunnel, and the calculation results were compared with those of the exact solution (analytical solution) and the classical FEM, respectively. The uniqueness, objectivity and rationality of the application of the FEM with grey variables in geotechnical problems are demonstrated.

Grey property of rock and soil system: A non-essential grey system
Using mathematical models to describe the rock and soil system is challenging, because mathematical models need detailed original information. However, the actual rock and soil system is complex, and full access to its information is limited based on current engineering technology and equipment. Therefore, the rock and soil system is a non-essential grey system. The grey property of rock and soil system is mainly reflected in the following aspects:

Complexity of rock and soil mass
Because the rock and soil system is a natural geological body meaning that its spatial distribution varies with crustal tectonic movement and metamorphism, the rock and soil system is heterogeneous. The stress field of rock and soil includes gravity stress, tectonic stress, seepage stress and temperature stress. There are many factors affecting the stress field, such as lithology, faults, groundwater, temperature, and etc. These factors are extremely irregular due to the heterogeneity and randomness of geological bodies. Therefore, the stress field is complex, and the boundary conditions of the studied area are difficult to determine.

The limited access to raw data
Due to the limitations of manpower, time, cost and field conditions, only partial geotechnical information can be obtained. The structure of rock and soil system is very complex, and most of them are in the blind area of people's sight. The data of geotechnical system collected by various engineering exploration methods are always limited and partial. The understanding of this actual system is incomplete, and it makes the rock and soil system a non-essential grey system.

Definition of grey number
Grey number (GN) is the most fundamental concept of grey system, a basic unit of grey system and grey mathematics.
Definition 3.1 [23]. Define P(c) as a connotation-cognition problem that can be quantified and I(c) as the cognition-information field of P(c). R is the field of real number and D is a subset of R. Suppose that � is an uncertain number about P(c), d o is the only potential true number of the uncertain number � and the acquiescing number�, and � satisfies the following condition,

Operation properties of grey number
The operation between grey numbers is usually the operation between number-covered set when engineering information is inadequate. Definition 3.2 [23]. Suppose that � i and � j are two continuous grey numbers (CGN), and their number-covered sets are Then we have that D ij is the number-covered set of � ij and it is as below: is the chaos of �, which represents the uncertainty degree of grey number [24], and is further used to evaluate the quality of the proposed grey models in Section 6.2.
Theorem 3.1 [23]. Supposing that � is a GN and D is its number-covered set. Let � 2 {+,−,×, �} be a binary operation. If a is a real number, then a� is GN, and its number-covered set is as below: Where 0 = 2 D j when the operator is "�" [24].

Definitions of grey parameters and variables, establishment of basic equations and differential equations with grey variables
Elastic mechanics is an important branch of solid mechanics. It studies the deformation and internal force of elastic objects under the action of external forces and other external factors, and it is also the basis of geotechnical problems. As studying the distributed parameter models, mechanical parameters, boundary conditions and stress conditions are essential.

Definitions of grey parameters and variables
The following analysis started from greying out basic variables and parameters, such as σ (stress), ε (strain), u (displacement), E (elastic modulus), μ (Poisson's ratio), because they are the foundation of basic elastic mechanics equations [25].
, where s ij and s ij represent the lower and upper bounds of grey stress at any point in space, respectively, i, j = x, y, z.

Establishment of differential equations with grey variables
The established differential equations with grey variables below satisfy the basic assumptions of elastic mechanics (continuity, uniformity, isotropy and small deformation) [25], and are used to study the equilibrium and deformation of any element of the elastic body.

The equilibrium equations with grey variables.
The equilibrium equations below with grey variables for two-dimension and three-dimension system were derived [25]. As shown in Fig 1, the grey load �F was decomposed in the x, y, z direction as (�f x �F y �F z ), according to the force equilibrium: Equilibrium equations with grey variables in x-axis direction can be obtained from Eq (1):

PLOS ONE
Finite element method (FEM) with grey system theory Similarly, equilibrium equations with grey variables in y-axis and z-axis direction can be obtained, respectively: @�s xy @x þ @�s yy @y þ @�s zy @z þ �f y ¼ 0 ð3Þ For a two-dimension system where the force and deformation of the object in the z-axis direction are ignored, the form of the equilibrium equations with grey variables can be further simplified as: @�s xy @x þ @�s yy @y þ �f y ¼ 0 ð6Þ

The geometric equations with grey variables.
The geometric equations with grey variables are discussed: The displacement components of the deformable body in the x-axis and y-axis directions are represented by �u(x, y, z) and �v(x, y, z). By using Taylor's formula, and ignoring the second-order and above terms and the deformation in the z-axis direction (Fig 2), equations below can be obtained: According to the definition of strain,�ε , small-deformation theory, and Eq (7), it can be obtained: Similarly, it can be obtained: According to the small-deformation theory, as the deformation angle is small, �γ x � tan �γ.

PLOS ONE
Finite element method (FEM) with grey system theory Therefore, the grey shear strain is: To sum up, the Eqs (11), (12) and (16) are the geometric equations with grey variables of the two-dimension problems. They can be written in matrix form as:

The physics equations with grey variables.
The physics equations with grey variables for two-dimension problems can be obtained from the stress-strain relationship determined by experiments: Where �τ xy is grey shear stress. They can be written in matrix form as: where �D ¼ �E is elastic matrix with grey variables.

Differential equations with grey variables.
By introducing the above geometric equations with grey variables into the physical equations with grey variables, the grey stress

PLOS ONE
Finite element method (FEM) with grey system theory represented by grey displacement can be obtained as: �s yy ¼ �E 1 À �m 2 ð @�v @y þ �m @�u @x Þ ð23Þ Substituting the above three Eqs (22)(23)(24) into the equilibrium Eqs (5 and 6) with grey variables, analytic equations for displacement can be obtained: With known boundary conditions, the Eqs (25) and (26) can be solved, and the analytical solutions of the grey displacements in the x-direction and y-direction can be calculated. In theory, the grey strain and stress can be obtained through the geometric equations with grey variables and physics equations with grey variables, respectively. But in real engineering cases, because of the complexity of boundary conditions, the equations are challenging to solve. In contrast, discrete numerical analysis is often more efficient. In the following discussion, the FEM is used to study the solution of the differential equations with grey variables. The core content is to convert the continuous infinite degree of freedom problem into a discrete finite degree of freedom problem for solving [26].

Finite element numerical solution method for differential equations with grey variables
The procedure of conducting the analysis of FEM with grey variables includes: firstly, the studied area was discretized into small elements; and then with the boundary conditions, load conditions and material properties with grey variables, the element stiffness matrix with grey variables, the assembled global stiffness matrix with grey variables and the global equation with grey variables were established. In the above process, the operation was completely carried out according to the operation properties of grey number. Finally, the results of grey displacement, grey stress and grey strain can be calculated. Among them, the establishment of element stiffness matrix with grey variables is the core of the steps [27]. Taking a 4-node rectangular element as an example, the establishment of the element stiffness matrix with grey variables is discussed.

PLOS ONE
Finite element method (FEM) with grey system theory As shown in Fig 3, the shape function [27] of 4-node rectangular element is:

PLOS ONE
Finite element method (FEM) with grey system theory where r and s are the coordinates in the parent unit. The shape function is used to connect the grey displacement of any point in the element with that of the element node: where �δ e = (�u 1 �v 1 �u 2 �v 2 �u 3 �v 3 �u 4 �v 4 ) T is grey displacement of the element node, x and y are the coordinates of sub unit, that is, the coordinates in the actual element. It can be obtained from Eqs (17) and (27): �ε ¼ @�u @x @�v @y @�v @x þ @�u @y where B ¼ According to the theory of FEM and virtual work principle, element stiffness matrix with grey variables can be derived by: where J ¼ @x @r @y @r @x @s @y @s The element stiffness matrix with grey variables can be then obtained by solving the Eq (29) through Gaussian Integral.
It can be found from the above process that geometric matrix and Jacobian Matrix are variables related to the shape function (N) and coordinates (x i , y i ). The grey variables only exist in grey elastic matrix. According to grey mathematics [28], the elastic matrix with grey parameters can be calculated by the operation of grey number. The stiffness matrix of each element with grey variables is further calculated. Finally, the global stiffness matrix with grey variables is assembled by the coding method: According to the theory of FEM [26], the global stiffness equation is: Where �δ = (�u 1 �u 2 . . . �u n ) T is global node grey displacement array, �u i = (�u i , �v i ) T represents the grey displacement of node i in two coordinate directions. �P = (�P 1 �P 2 . . . �P n ) T is global node grey load array, �P i = (X i Y i ) T represents the grey load of node i in two coordinate directions. From one-to-one correspondence, it can be seen from the relationship between the variables that the physical meaning of global stiffness matrix with grey variables is the force at the i-th node of the element caused by the unit deformation of the j-th degree of freedom of the element. The above process was programmed by Matlab. The calculation of FEM with grey variables was carried out according to the operation properties of grey number.

Analytical solution and classical FEM solution
A deep-buried circular mining tunnel was selected as an example to study the numerical method with grey variables, where a stability evaluation of the surrounding rock is required. The surrounding rock of the tunnel is composed of sandstone and shale, and the geological structure is very complex, with heterogeneity and anisotropy. This geotechnical problem has many grey system characteristics, therefore, the FEM with grey variables is used to solve the problem. For clearly study, a section is chosen and it is assumed to be elastic. The deep circular tunnel model with uniform stress field of 100m � 100m was built and shown in Fig 4. According to the classical FEM, the parameters and variables are treated to be deterministic, meaning that they have certain numerical values. The parameters in Table 1, are mean values of the measured or experimental ones in the example. However, the determined certain values can hardly cover all the information. In the example, the elastic modulus (E) of surrounding rock mass is 1 � 10 8 N/ m 2 , and the Poisson's ratio (μ) is 0.3. The tunnel's radius is 3m, and the surrounding area bears uniform initial stress, which is 2 � 10 6 N/m 2 . Supporting force is applied to the inside of the tunnel, the direction is perpendicular to the tangent direction, and the magnitude is 1.06 � 10 6 N/m 2 . The displacement at the boundary of 100m was 0.0011m as calculated [25].
The pre-processing was completed by FEM software Abaqus, which included the establishment of the model, the input of loads and boundary conditions, and the division of the meshes. The divided grid was shown in Fig 5 with 227 elements and 255 nodes in total. The output data from Abaqus (including the node number, node coordinates and load of each node of each element) were used as the input information for the FEM analysis with the grey variables.
When dividing the grid, as the size of the mesh becomes more refined, the numerical solution will approach the exact solution (analytical solution), meaning the numerical solution is getting more accurate. However, small mesh sizes will lead to an excessive increase in the computational cost. Therefore in this study, refined mesh size was only adopted in the key areas in order to meet the accuracy, while for the non-key areas (such as the center distance greater than 10 m in this example), the grid division was relatively rough in order to improve the computational efficiency, as shown in Fig 5. In order to ensure the correctness of the FEM with grey variables, the classical FEM method was first programmed, and the data pre-processed by Abaqus were used as its input information for calculation. The results were compared with the analytical solution results. Ten representative nodes with different center distances were selected, as shown in Table 2. The comparison of displacement solutions of the nodes was shown in Fig 6. It can be seen from Table 2 that the mesh size increases with the increase of center distance, while the accuracy of

PLOS ONE
Finite element method (FEM) with grey system theory numerical solution gradually decreases (i.e. the relative error increases). Regarding the region that we focus on (center distance 3-5 m, where the deformation is critical to tunnel's support and operation), the accuracy of numerical solution is satisfying, and the relative error is about 1% (i.e. node 69 and node 160). The calculated displacement by the classical FEM for the meshed grid above was plotted in the centrally symmetric displacement contour maps (Figs 7 and 8). It can be seen from Fig 7 that as the center distance increases, the change of displacement becomes smaller. The displacement of the center distance range of 10m-20m is 0.005m, and the displacement of the center distance ranging from 20m to 100m is about 0m. Therefore, we will focus on studying

PLOS ONE
Finite element method (FEM) with grey system theory the displacement within the center distance of 10m, as shown in Fig 8. The related stress and strain were not discussed in this paper.

Determination of optimal solution of FEM with grey variables and its comparison with the analytical solution and classical FEM numerical solution
The determination of parameters and variables in the above example are changed in sense of that the max value and min value of parameters and variables make up the number-covered

PLOS ONE
Finite element method (FEM) with grey system theory set of each grey parameters and grey variables, and in this way, all information will be included.
is the grey displacement to solve. The calculation process is computationally intensive, and the results will be distorted if only following the solution method of the interval grey number equations (Fig 9). In this regard, a study on the solution of the grey linear equations was conducted. It can be seen from Eq (30) that �δ is inversely proportional to �K and proportional to �P. The number-covered sets of �δ can be calculated by the following equation: The Eq (29) shows that when the element stiffness matrix with grey variables is calculated, only the elastic matrix �D contains grey parameters. When �D is calculated, according to the operation properties of grey number, the lower bound of each element obtained is less than the upper bound. However, when the element stiffness matrix with grey variables is calculated, the lower bound of some elements becomes greater than the upper bound. The load matrix with grey variables also has the same problem. By rearranging the elements of the stiffness matrix with grey variables and the load matrix with grey variables with the rule that the lower bound is smaller than the upper bound, and calculating with Eq (31), the results can be derived as the displacement contour map in

PLOS ONE
Finite element method (FEM) with grey system theory It can be seen from the comparison in Fig 9 that most of the displacement of Umax is greater than that of Umin, but the displacement distribution of Umax is irregular making its reference value limited. Although the displacement distribution of Umin is more regular, its reference value is also poor because it is contrast to the above result that the displacement is inversely proportional to center distance. This is because after the element stiffness matrix with grey variables and the load matrix with grey variables were rearranged, the originally uniform and regular solution area became irregular. Therefore in this study, after calculating the elastic matrix �D (in which the lower bound of each element in �D is less than the upper bound) through the operation properties of grey number, and further calculating the element stiffness matrix with grey variables (some lower bounds in the stiffness matrix of the grey element may be greater than the relative upper bounds), the element stiffness matrix with grey variables and the load matrix with grey variables were not rearranged. The obtained displacement solutions were plotted in Fig 10. It can be seen from Fig 10 that all the displacement of Umax is greater than Umin, and the grey displacement distribution is regular. Comparing with the classical FEM, it can be found that the displacement results of the classical FEM are included within the chaos (defined in Definition 3.2 (d)) of grey displacements of the FEM with grey variables, which implies the model from Eq (31) is more reasonable and has more information considered than the classical FEM. For example, when the center distance is 4m, the displacement of the classical FEM is 0.025m, and the displacement of FEM with grey variables is [0.02, 0.04] m.
Similarly, the other solution can be calculated from Eq (32) (in Fig 11) when the elastic matrix �D (in which the lower bound of each element in �D is less than the upper bound) is calculated through the operation properties of grey number, and the element stiffness matrix with grey variables and the load matrix with grey variables are not rearranged: It can be seen from Fig 11 that the Umax is greater than the Umin, the distribution is regular, and the displacement results of the classical FEM are also included within the chaos of grey displacements of the FEM with grey variables, making the results calculated by Eq (32) correct and reasonable. Furthermore, the chaos of grey displacements calculated by Eq (32) are included within those calculated by Eq (31). Therefore, Eq (31) is more advantageous than Eq (32), and it provides the optimal solution of FEM with grey variables in this study.
The displacement solution (chaos of grey displacement) of FEM with grey variables was compared with the analytical displacement solution and the classical FEM displacement solution, respectively. Ten nodes were selected as mentioned above, and the displacement comparison was plotted in Fig 12. It can be seen that the analytical displacement solution and the classical FEM displacement solution are both included in the chaos of grey displacement solution of FEM with grey variables, indicating that the FEM with grey variables is reasonable. Moreover, bigger chaos of grey displacement is obtained when center distance of the node is smaller (meaning a greater grid accuracy, i.e. node 69 and node 160). This is because smaller center distance corresponds to bigger stress, and further causes more complex deformation situation and larger chaos of grey displacement. It is worth pointing out that greater chaos of grey displacement does not represent greater calculation error, on the contrary it means the results are more safe by including more uncertainty information which is missed in the classical FEM.

1.
A new idea of numerical simulation of geotechnical problems is proposed, i.e. constructing the numerical model with grey system theory. Based on the grey system theory, this paper https://doi.org/10.1371/journal.pone.0270400.g010

PLOS ONE
Finite element method (FEM) with grey system theory analyzed the grey properties, grey parameters and grey variables of rock and soil system, established the basic equilibrium equations and the FEM equations with grey variables, and used the operation properties of grey number to solve the FEM equations with grey variables. https://doi.org/10.1371/journal.pone.0270400.g011

PLOS ONE
Finite element method (FEM) with grey system theory 2. An example of a deep-buried circular mining tunnel was executed, and the displacement solution of FEM with grey variables was obtained and compared with the analytical displacement solution and classical FEM displacement solution. The comparison results show that both the analytical displacement solution and classical FEM displacement solution were included in the chaos of grey displacement calculated by FEM with grey variables, and the accuracy of calculation results was improved with more refined mesh size. This demonstrates that the results of the proposed grey numerical model contain more information (i.e. the bigger chaos of the displacement), and the potential danger of losing some information

PLOS ONE
Finite element method (FEM) with grey system theory in case of applying the classical numerical model can be avoided. So the proposed grey numerical model is more reliable when dealing with real geotechnical problems.
3. The chaos of grey displacement of FEM with grey variables increases as the center distance decreases, due to the fact that the smaller center distance of the node (meaning higher the grid accuracy, i.e. node 69 and node 160) is, the bigger stress and the more complex deformation situation is. A larger chaos of grey displacement suggests that the result be more reliable, making the proposed model capable of handling the complexity in real cases.